Effect of helicity and rotation on the free decay of turbulent flows 
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The self-similar decay of energy in a turbulent flow is studied in direct numerical simulations with 
and without rotation. Two initial conditions are considered: one non- helical (mirror-symmetric), 
and one with maximal helicity. The results show that, while in the absence of rotation the energy in 
the helical and non-helical cases decays with the same rate, in rotating flows the helicity content has 
a major impact on the decay rate. These differences are associated with differences in the energy 
and helicity cascades when rotation is present. Properties of the structures that arise in the flow at 
late times in each time are also discussed. 



Turbulence is ubiquitous in nature, and many turbu- 
lent flows are also rotating. The effect of rotation be- 
comes important when the Rossby number (the ratio of 
the convective to the Coriolis acceleration) is sufficiently 
small. Mid-latitude synoptic scales in the atmosphere, 
stellar convective regions, and turbomachinery are ex- 
amples of such flows. Helicity (alignment of the velocity 
and the vorticity) is also important for many processes 
in astrophysical, geophysical and engineering flows. As 
an example, helical flows were proposed as the reason for 
the stability of rotating convective thunderstorms [l[ ■ 

Studies of isotropic and homogeneous helical turbu- 
lence 0, showed that both the helicity and the en- 
ergy are transferred toward smaller scales with constant 
fluxes. Moreover, it was observed that the scaling of 
the energy (Kolmogorov's law) was unchanged by the 
presence of helicity. As a result, helicity is expected to 
globally arrest the energy decay, but not to change its 
self-similar decay rate. Non-helical rotating turbulence 
has been studied in detail (see e.g., [|| and references 
therein) , but perhaps because of the similarities between 
the helical and non-helical cases in non- rotating flows, 
not much attention has been paid to helical rotating tur- 
bulence. 

The lack of detailed studies of rotating helical flows is 
remarkable considering the relevance of both helicity and 
rotation in many astrophysical and geophysical processes. 
In this work, we study the effect of rotation and of helic- 
ity in the self-similar decay of energy in turbulent flows. 
Even in isotropic and homogeneous turbulence, the law 
for the decay rate of energy is a matter of debate [J] . It 
is known that it depends on properties of the infrared en- 
ergy spectrum (i.e., the spectrum at scales larger than the 
energy containing scale), and may depend on other sta- 
tistical properties of the initial conditions. As a result, in 
this letter we will consider only two flows with the same 
infrarrcd spectrum and the same energy decay rate in the 
absence of rotation, and study how the presence of he- 
licity and rotation changes their decay. New decay laws 
are found for helical rotating flows in numerical simula- 
tions. The results are then interpreted in terms of how 
the energy and helicity cascades are modified when rota- 
tion is present, and a phcnomcnological theory that is in 



TABLE I: Parameters used in the simulations: v is the kine- 
matic viscosity, Q is the rotation rate, Re is the Reynolds 
number, Ro is the Rossby number, Ro" is the micro-Rossby 
number, Ek is the Ekman number, and h is the relative helic- 
ity of the initial conditions. The values of Re, Ro, Ro" , and 
Ek are given at the time of the peak of dissipation t* . 
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agreement with the simulations is finally discussed. 

The numerical simulations solve the Navier-Stokcs 
equations for an incompressible fluid in a rotating frame, 
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where u is the velocity field (V-u = 0),u; = Vxu 
is the vorticity, V is the total pressure (modified by the 
centrifugal term), and v is the kinematic viscosity. We 
chose the rotation axis to be in the z direction, $1 = Oz, 
with Q the rotation frequency. Our integration domain is 
a periodic box of length 2tt. Two sets of runs were done 
at resolutions of 256 3 (set A) and 512 3 grid points (set 
B) using a pseudo spectral code. The parameters for all 
the runs arc listed in Table |H 

To simulate systems with different amount of relative 
initial helicity (h = H/ (|u||u;|) where H = (u • u>) is 
the flow net helicity), two flows were considered as ini- 
tial conditions: the Taylor-Green (TG) flow and the 
Arn'old-Beltrami-Childress (ABC) flow @. The TG flow 
is non-helical, and has zero energy in the k z = mode, 
whose amplification observed in the rotating cases (see 
below) is thus only due to a cascade process. The TG 
flow was chosen for its importance in hydrodynamics; it 
was originally motivated as an initial condition which 
leads to rapid development of small spatial scales. It 
also mimics the von Karman flows between two counter- 
rotating disks used in several experiments. The ABC 
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flow is an eigenfunction of the curl operator and as a 
result has maximum helicity. It was used as a paradig- 
matic example to study helical flows in the atmosphere 
[!|. Both flows develop, after a short time, an infrared 
energy spectrum proportional to k 2 ; this is important to 
ensure we can compare the decay rates of both flows. 

The simulations were started using a superposition of 
these flows from wavenumbers k = 4 to 14. Runs with 
zero relative helicity have a supperposition of TG flows 
and runs with h w 0.95 have a supperposition of ABC 
flows. The initial energy spectrum from k — 4 to 14 
was proportional to fc -4 , leaving enough spectral space 
to allow for direct and inverse transfer of energy. All the 
runs were extended for over 40 turnover times, and the 
dissipative range was properly resolved. Times in Table 
Uand in the figures are expressed in units of the turnover 
time at t = 0, T = L/U, where L « 27r/fc = 27r/4 is the 
initial integral scale of the flow and {/ ~ 1 is the initial 
r.m.s. velocity. 

Several Reynolds, Rossby, and Ekman numbers can be 
defined for the runs. We will consider here the Reynolds 
number based on the integral scale Re = UL/v, and 
the accompanying Rossby number Ro — U/(2CIL). The 
integral scale is defined as L = 2ir J E(k)k~ 1 dk/E where 
E(k) is the isotropic energy spectrum and E is the total 
energy. The Ekman number is then defined as E^ = 
Ro/Re. At t = runs in set A have Re = 1047 and runs 
in set B have Re = 2244. At the same time the Rossby 
number is Ro = 0.08 for runs A3, A4, and B1-B2; the 
initial Ekman number is then E k = 7.6 x 10~ 5 in runs 
A3 and A4, and E k = 3.6 x 10~ 5 in runs B1-B2. Values 
of the controlling parameters for each run at the time of 
maximum dissipation t* are given in Table |TJ 

It is also convenient to introduce a micro-scale Rossby 
number as Ro^ = lu/(2Q) Q- This number can be in- 
terpreted as the ratio of the convective to the Coriolis 
acceleration at the Taylor scale, a scale characteristic of 
the turbulent inertial range. In order for anisotropics to 
develop in the simulations, the Rossby number Ro must 
be small enough for the rotation to affect the free decay- 
ing turbulence, but the micro- Rossby Ro u must be larger 
than one for scrambling effects of inertial waves not to 
completely damp the nonlinear terms (which leads to a 
pure exponencial viscous decay) p|. We are thought in- 
terested in simulations with moderate rotation rates to 
ensure Ro < 1 and with Re large enough to have enough 
scale separation between the energy-containing and the 
dissipative scale. Note in Table fl] how Ro and Ro u arc 
one order of magnitude apart in all runs at time t* . This 
difference is sustained in time in the runs during the self- 
similar energy decay. 

Figure [T] shows the time history of the energy and cn- 
strophy in runs Al, A2, B1-B2. In all runs there is a 
self-similar decay after the time t*. Runs A3 and A4 
(not shown) decay respectively as runs Bl and B2, al- 
though the time span of the self-similar stage is shorter. 
In both runs with O = 0, the energy decays as ~ t~ 2 . As 
will be discussed later, this exponent corresponds to the 
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FIG. 1: Energy evolution for runs Al (solid), A2 (dotted), Bl 
(dashed) and B2 (dot-dashed). Runs Al and A2 have SI = 
and decay with the same rate independently of the helicity 
content. When rotation is turned on the decay rates differ. 
Inset: enstrophy evolution for runs Al and A2. 



decay of a flow with constant integral length [9j]. 

The presence of helicity in run A2 doesn't seem to 
affect the self-similar decay of turbulence, as predicted in 
Ref. [ltj. This is in good agreement with the fact that 
helicity doesn't change the spectral index of the energy 
in non-rotating turbulence [2J. However, the self-similar 
decay in run A2 starts at a later time, as also previously 
reported in Ref. [Il|. This is associated with the slow- 
down in the generation of small-scales in helical flows [ll] , 
which results in a longer time to reach the maximum of 
dissipation in run A2 (see the inset in Fig. [T]). 

In runs with rotation (A3, A4, Bl and B2), a transient 
is also observed before t*. Then, self-similar decays with 
different power laws for the energy are found in all runs. 
Runs A3 and Bl have a decay near i -1 [H|,[l4[, while the 
runs with maximum helicity (A4 and B2) follow a decay 
slightly faster than t -1 / 3 . Contrary to the non-rotating 
case, we can appreciate now how the presence of initial 
helicity changes the self-similar energy decay. 

The time evolution is accompanied by a change in the 
shape of the energy spectrum (see Fig. EJa)). While 
the runs without rotation show at t* and during the self- 
similar decay an spectrum consistent with Kolmogorov 
scaling in both the energy and the helicity, runs with ro- 
tation show different scaling laws. Run Bl (which has a 
larger scale separation than run A3) shows an isotropic 
energy spectrum compatible with k (the anisotropic 
energy spectrum is also compatible with a kj 2 scaling) 
[13, Ell]- On the other hand, the energy spectrum in run 
B2 is slightly steeper than in run Bl, while the helicity 
spectrum in this run is shallower than the energy spec- 
trum. The product of both spectra is consistent with a 
k~ A scaling law, as predicted using phenomenology for 
rotating helical flows in (l8j . 

In the runs with rotation, a change in the small-scale 
spectrum is observed at t w 20. As energy piles up at the 
largest available scale in the box and column-like struc- 
tures form in the velocity, the small-scale energy spec- 
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FIG. 3: Ratio of the energy in modes with fen = to the 
total energy in runs Bl (dashed) and B2 (solid). Despite the 
different initial conditions, the curves grow monotonically to 
f, showing a transfer of energy to perpendicular modes. Inset: 
time evolution of the Shebalin angle for the same runs. 



FIG. 2: (a) Isotropic energy spectrum for run B2 at times 
t « 2 (thin), t « 12 (thick) and t ~ 30 (dot-dashed). At 
early times the spectrum follows a ~ fc~ 2 ' 5 power-law (slope 
shown as a reference) and becames later dominated by the 
piling up of energy in the large scales. Inset: Isotropic energy 
spectrum for run B2. A k~ 2 slope is plotted for reference, (b) 
Energy (thin) and helicity (thick) normalized fluxes for run 
B2 at times t « 5 (solid) and t ~ 13 (dashed) respectively. 
Inset: Energy flux for run Bl at times t ~ 2 (solid) and t m 11 
(dashed). 



trum becomes steeper. This is accompanied by a decrease 
in the decay rate of the energy (see Fig. [IJ . This process 
is reminiscent of the change observed in the free-decay 
of two-dimensional turbulence, when the coalescence of 
large-scale vortices at late times leads to a steeper energy 
spectrum and a change in the self-similar decay [IB] 

Runs Al and A2 developt direct energy fluxes toward 
small-scales while the rotating non-helical runs (A3 and 
Bl) show both a direct and an inverse energy cascade 
(see the inset of Fig. [2] (b)). In the helical runs (Fig. [2] 
(b)), at t w 2 we observe maximum flux of energy and 
helicity toward smaller scales, evidenciating both energy 
and helicity have a direct cascade. However, the helicity 
flux is larger than the direct energy flux. Later, an in- 
verse cascade of energy can be clearly identified from the 
negative energy flux at large scales. At t w 13, the co- 
existence of both an inverse cascade of energy and direct 
cascades of energy and helicity is observed. 

Although Figs. [2] (a) and (b) show the isotropic spectra 
and fluxes, the flows in the rotating case are anisotropic. 
Most of the energy in the spectra is in modes perpen- 
dicular to the axis of rotation, and the anisotropic spec- 
turn E(k±) and flux II(fcjJ look similar to the ones pre- 
viously discussed. We present instead some global in- 
dications of the development of anisotropies. Figure [3] 
shows the ratio of the energy in perpendicular modes 



(i.e., modes with k\\ = 0) to the total energy in runs 
Bl and B2. Differences at early times are due to diffcr- 
emt initial conditions. However, it can be seen that both 
curves grow monotonically to a value near 1, indicating 
the flows evolve toward anisotropic states as the energy 
is transferred to perpendicular modes Q ■ 

A measure of anisotropies in the small-scale fluctua- 
tions is given by the Shebalin angle, defined as 

tan2g = 2 Efer^) ; (2) 

The evolution of this angle in runs Bl and B2 can be seen 
in the inset of Fig. [3J It grows monotonically although in 
run B2 it reaches a maximum at t « 12 and then seems 
to saturate. A similar behavior is observed for L|| (not 
shown) which again grows monotically in Bl but reaches 
a maximum in B2 at aproximately the same time. In 
all cases, tan 2 6* >> 2 which corresponds to anisotropic 
flows. 

The increase of the correlation lengths, together with 
the growth of E(ku = 0)/E and of tan 2 8 speaks of 
anisotropization of the flows This tendency towards two- 
dimensionalization is confirmed by exploring the flows in 
real space. Fig. U shows visualizations of the r.m.s. vor- 
ticity with superimposed velocity field lines for runs A3 
and A4 at late times (t 45). In both cases, a strong 
anisotropy is observed with large scale column-like struc- 
tures in the vorticity (similar structures are observed if 
energy density is visualized instead). However, a signifi- 
cant difference in the geometry of the flow in the columns 
is observed between runs A3 and A4 . While the flow in 
the columns of the helical run is strongly helical (the flow 
goes either up or down in the entire column, thus giving 
helical velocity field lines), the columns of the non-helical 
run have no net helicity (the flow goes up and down in- 
side the same column). This indicates that even at late 
times, the properties of the emerging structures in rotat- 
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FIG. 4: (Color online) Visualizations of the r.m.s. vorticity 
for t m 45 with superimposed velocity field lines for runs A3 
(right) and A4 (left). Cuts of v z and H on a plane perpen- 
dicular to ft are shown. The color table gives the amplitudes 
normalized to the maximum of each quantity. 

ing flows depend on the initial helicity content. 

The results presented so far show a novel and distinct 
evolution in the free decay of helical and non-helical flows 
when rotation is present. These differences can be un- 
derstood in terms of a simple phenomenological theory. 
For the non-rotating cases we make use of the classi- 
cal Kolmogorov phenomenology which leads to the well 
known energy spectrum E(k) ~ e 2 / 3 fc -5 / 3 , where e repre- 
sents the energy dissipation rate. From the Navier-Stokes 
equation it is easy to show that for a freely-decaying flow 
dE/dt ~ E 3 / 2 / L which leads to a self-similar decay law 
E(t) ~ t~ 2 if the integral scale is assumed constant 0,01. 

In the case of rotating non-helical flows, it is often as- 
sumed that E(k) ~ e 1 / 2 ^ 1 / 2 /^ 2 (for simplicity we use 



the isotropic wavenumber k, although the arguments here 
can be easily generalized to the anisotropic case replac- 
ing A; by k±). This spectrum was observed in simula- 
tions [TtJ , and obtained from closures [l4| . It can also be 
derived from phenomenological arguments assuming the 
inertial waves slow down the energy cascade [l9j . Using 
this spectrum it follows that dE/dt ~ (E/L) 2 resulting 
in a decay law E(t) ~ t -1 . 

The novel case of rotating helical flows differs from the 
previous two cases. The direct transfer is dominated by 
the helicity cascade. Writting the helicity flux as S ~ 
hi/(£lrf) where hi is the helicity at the scale I and t; the 
eddic turnover time, leads to spectra E[k) ~ k~ n and 
H(k) ~ k n ~ 4 , where n = 5/2 for the case of maximum 
helicity [l8||. In this case, dimensional analysis leads to 
E{t) ~ e 1 / 4 ^ 5 /^" 5 / 2 and to E{t) ~ t' 1 / 3 . Note that for 
flows with initial helicity between zero and the maximum, 
the decay rate is between —2 and —1/3. A run with 
relative helicity of h m 0.4 was done to verify this. 

The results indicate that helical and non-helical rotat- 
ing flows differ in their scaling laws and decay rates. This 
is different than the case of non-rotating flows where he- 
licity doesn't change the energy scaling. The observed 
differences indicate that helical rotating flows must be 
studied in detail, and open new directions in the study of 
rotating turbulence. This can also be the starting point 
to elucidate the actual role helicity plays in the decay of 
flows in nature where rotation is often relevant 
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